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Abstract 

We study the 2-dimensional Navier-Stokes-Poisson equations with density-dependent vis- 
cosity 9 = 1/2 without pressure of gaseous stars in astrophysics. The analytical solutions with 
collapsing in radial symmetry, are constructed in this paper. 

1 Introduction 

The evolution of a self-gravitating fluid can be formulated by the Navier-Stokes-Poisson equations 
of the following form: 

{Pt+V. (pu) =0, 
(pu)j-l-V • (pu 'Siu) + VP ==-pV$ -I- vis{p, u), (1) 

where a (A'') is a constant related to the unit ball in i?^: a{l) = 2; a (2) = 27r and For > 3, 

aiN) ^ NiN - 2)ViN) ^ N{N - 2) ^^^^^ ^ , (2) 

where V{N) is the volume of the unit ball in and F is a Gamma function. And as usual, 
p — p{t,x) and u = u{t,x) E are the density, the velocity respectively. P = P{p) is the 
pressure. 
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In the above system, the self-gravitational potential field $ = x) is determined by the 
density p through the Poisson equation. 
And vis{p,u) is the viscosity function: 

vis(p,u) = Sj(p{p)sj»u). (3) 

Here we under a common assumption for: 

pip) - KpO (4) 

and K and 9 > are the constants. In particular, when = 0, it returns the expression for the u 
dependent only viscosity function: 

vis{p,u) = kAu. (5) 

And the vector Laplacian in u{t, r) can be expressed: 

N ~l N -I 

Au = Urr H Ur j U. (6) 

The equations and ^2 {vis{p,u) 7^ 0) are the compressible Navier-Stokes equations with 
forcing term. The equation ([Ija is the Poisson equation through which the gravitational potential 
is determined by the density distribution of the density itself. Thus, we call the system ([ij the 
Navier-Stokes-Poisson equations. 

Here, if the vis{p,u) — 0, the system is called the Euler-Poisson equations. In this case, the 
equations can be viewed as a prefect gas model. For = 3, ([T]) is a classical (nonrelativistic) 
description of a galaxy, in astrophysics. See [5], [5] for a detail about the system. 

P = P{p) is the pressure. The 7-law can be applied on the pressure P{p), i.e. 

P{p) = Kp'^ = ^, (7) 
7 

which is a commonly the hypothesis. The constant 7 = cp/cy > 1, where cp, are the specific 
heats per unit mass under constant pressure and constant volume respectively, is the ratio of the 
specific heats, that is, the adiabatic exponent in ([71)- In particular, the fluid is called isothermal if 
7 = 1. With K = 0, we call the system is pressureless. 

For the 3-dimensional case, we are interested in the hydrostatic equilibrium specified by u = 0. 
According to 2J, the ratio between the core density p{0) and the mean density pfor6/5<7<2 
is given by 

^ ^'%{^)] (8) 



p(o) V ^ 



z=zn 



where y is the solution of the Lane-Emden equation with ^=1/(7— 1), 

y{z) + -y{z) + = 0, y(0) = a > 0, y(0) = 0, n = (9) 

z 7 ^ -L 
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and zq is the first zero of y(zo) = 0. We can solve the Lane-Emden equation analytically for 



Vanaliz) = ' 



l-lz^, n = 0; 
sinz 

' 1 

-., n — 0, 



(10) 



and for the other values, only numerical values can be obtained. It can be shown that for n < 5, 
the radius of polytropic models is finite; for n > 5, the radius is infinite. 

Gambin [3] and Bezard [Ij obtained the existence results about the explicitly stationary solution 
{u — 0) for 7 = 6/5 in Euler-Poisson equations: 



P-['^] + (11) 



where A is constant. 

The Poisson equation ^3 can be solved as 



$(i,a;)=/ Gix-y)p{t,y)dy, (12) 
Jr." 

where G is the Green's function for the Poisson equation in the A^-dimensional spaces defined by 

c|, iV = l; 

G{x) = { log|x|, iV = 2; (13) 

In the following, we always seek solutions in radial symmetry. Thus, the Poisson equation ((ija is 
transformed to 

r^-i$„(t,a;) + {N - 1) r^'Hr^a (N) pr'^-\ (14) 

In this paper, we concern the analytical solutions with core collapsing for the 2-dimcnsional 
pressureless Navier-Stokes-Poisson equations with the density-dependent viscosity. And our aim is 
to construct a family of such core collapsing solutions. 

Historically in astrophysics, Goldreich and Weber [5] constructed the analytical collapsing 
(blowup) solution of the 3-dimcnsional Eulcr-Poisson equations for 7 = 4/3 for the non-rotating 
gas spheres. After that, Makino |7j obtained the rigorously mathematical proof of the existence of 
such kind of collapsing solutions. And in [3J, the extension of the above collapsing solutions to the 
higher dimensional cases {N > 4). After that, for the 2-dimensional case, Yuen constructed the 
analytical collapsing solutions for 7 = 1, in [H]. 

For the construction of the analytical solutions to the Navier-Stokes equations in , the 
Navier-Stokes-Poisson equations in without pressure with = 1 and in without pressure 
with 9 = 5/4, readers may refer Yuen's recent results in [9], [10], [11] respectively. 
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In this article, the analytical collapsing solutions are constructed in the pressureless Navier- 
Stokes-Poisson equations with density-dependent viscosity in R^, with 9 — 1/2, in radial symmetry: 
1 

Pt + Upr + pUr + -pu = 0, 

^ 2ttp r 11 (1^) 

p (Uf + UUr) ~ / p(t, s)sds + + {K.p^/'^){u„ H Uj. -u) , 

r Jq J. J.Z 

in the form of the following theorem. 

Theorem 1 For the 2- dimensional pressureless Navier-Stokes-Poisson equations with 9 — 1/2, in 
radial symmetry, f J5)) . there exists a family of solutions, 



y{z)+\y{^) - = 0' 2/(0) = " > 0, y(0) = 0, 



where T > 0, k > 0, C 7^ and a are constants. 

In particular, for C > 0, the solutions collapse in the finite time T/C . 



2 Separable Blowup Solutions 

Before presenting the proof of Thcorcm[l] we prepare some lemmas. First, we obtain the solutions 
for the continuity equation of mass in radial symmetry (|15p 1 . 



Lemma 2 For the equation of conservation of mass in radial symmetry: 



there exist solutions, 



Pt + Upr + pUr + -pu — 0, (17) 

r 



1 — C 

p(t,r) = 5 , uit,r) = —r, (18) 

^ ' {T-Ctfy{^Y ^' ' T-Ct ' ^ ^ 



with the form y ^ and y Q , C and T > are constants. 

Proof. We just plug ^ into HH). Then 

1 

Pt + Upr + pUr + -pu 

r 

{-2){-C) , (-2)(-l)(-C)rz/(y^) 



(T-Ct)%(y^)2 (r-Ci)%(T^)3 

{-C)r (-2) yij^t) , 1 hC) 



T - Ct [T - Ctfy{.j^^f T-Ct [T - Ctfyij^^Y T - Ct 
1 1 (-C) 
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The proof is completed. ■ 

Besides, we need the lemma for stating the property of the function y{z). The similar lemma 
was already given in Lemmas 9 and 10, ^8j, by the fixed point theorem. For the completeness, the 
proof is also presented here. 

Lemma 3 For the ordinary differential equation, 



1 

—1 

2/(0) = a > 0, 2/(0) = 0, 



2/(z)+-?/(z) " -rn = 



where a is a positive constant, 

has a solution y{z) E and lim y{z) — oo. 



Proof. By integrating ([T9)l . we have, 

1 

z Jq y{s) 

Thus, for < z < zq, y{x) has a uniform lower upper bound 



y{z) = - I 7732 sds > 0. (20) 



y{z) > 2/(0) = a > 0. 
As we obtained the local existence in Lemma [31 there are two possibilities: 

(1) 2/(2;) only exists in some finite interval [0, zq]: (la-) 1™ y{z) — 00; (lh)y{z) has an uniformly 

z—tzo- 

upper bound, i.e. y{z) < ao for some constant ao. 

(2) 2/(2;) exists in [0, +00): (2a) lim y(z) — 00; (2b)j/(z) has an uniformly upper bound, i.e. 

z — ' + 00 

y{z) < P for some positive constant /3. 

We claim that possibility (1) does not exist. We need to reject (lb) first: If the statement (lb) is 
true, becomes 

az a s , , ^ , 
ds > y{z). (21) 



2a2 z 

Thus, y{z) is bounded in [0, zq]. Therefore, we can use the fixed point theorem again to obtain a 
large domain of existence, such that [0, 2:0 + ^] for some positive number 6. There is a contradiction. 
Therefore, (lb) is rejected. 

Next, we do not accept (la) because of the following reason: It is impossible that lim y{z) — 00, 

z^zo^ 

as from (|21l) . y{z) has an upper bound in [0, zq]: 

S > y(^)- (22) 



Thus, (1221) becomes, 



rzo 

y{zo) = 2/(0) + / y{s)ds 
Jo 

fZn 

< Ct + -:-^ds 



a 



2a2 
2q;2 
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Since y{z) is bounded above in [0, zq], it contracts the statement (la), such that hm y{z) = oo. 
So, we can exclude the possibihty (1). 

We claim that the possibihty (2b) doesn't exist. It is because 



Then, we have. 



y(2)>a+^2'. 



(23) 



By letting z — > oo, (pS)) turns out to be. 



y(z) = oo. 



Since a contradiction is established, we exclude the possibility (2b). Thus, the equation exists 
in [0, +oo) and lim y(z) — oo. This completes the proof. ■ 

Here we are aheady to give the proof of Theorem [TJ 
Proof of Theorem 2. From Lemma [U it is clear for that (fTBl) satisfy (fl3|) i . For the momentum 
equation (fT5|) 9. we get, 



piut + UUr) ^ f pit, S)sds — \p{p)]rUr — p(p)(Urr + -Ur ttu) 

r J. j,z 



(_C)(-1)(-C) , (-C) (-C) 
■r + — —r 



(T - Ctf 



T-Ct T-Ct 



27rp 



1 



sds 



T-Ct 



■T-Ct r^T-Ct 



d_ 

dr 

1 
2 

2lTp 



(T-Ctfyi^^r 



1/2 



c 



1 -1/2 



(T-Ct)MT^)^ 
1 



27rp 
r 

-2) 



T-Ct r J {T-Ctfy{^J- 



-sds 



{T - CtYy{.j^^f T-Ct T-Ct 



r J {T-Ctfy{M^ 



-C pyijr^t) , 27rp 



sds 



^ {T- Ct) 

-P 



r J {T-Ctfy{M^ 



sds 



(T - Ct) 



(T - Ct) 



2tt 



K^T-Ct' r{T - Ct) J y{j^^)^ 

r/{T-Ct) 

-^ViTT; — 7^) / -^^TdT 



sds 



/i^^'T-Ct' 



T~Ct 



y{r)' 



(24) 
(25) 



(26) 



(27) 
(28) 
(29) 

(30) 
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By letting T^r/{T - Ct), it follows: 

-P 



[T - Ct) 
-P 



C V , 



27r 



TliT-Ct) 



T-Ct 



y{T)' 



-rdr 



Q 



{T-cty \T~Ct 

And denote z = r/(T — ct), 



(31) 



<3(5^) = 0(=)=;|.W-|/^ 



rdr. 







Differentiate Q{z) with respect to z, 



^^^^ ^ 7^^^^^ - ^ + ^7 



a y(z) 27r ? 1 



rR. z j ?;(t) 



1 



= — QW, 

z 

where the above result is due to the fact that we choose the following ordinary differential equation: 

y(0) = a > 0, y(0)= 0. 

With (5(0) = 0, this implies that Qiz) = 0. Thus, the momentum equation (fTS]) ? is satisfied. 
With Lemma |3] about ?;(z), we are able to show that the family of the solutions collapse in finite 
time T jC. This completes the proof. ■ 

The statement about the blowup rate will be immediately followed: 



Corollary 4 The collapsing rate of the solution \16\) is 



Jim p(t,0)(r-Cif> 0(f). 



(32) 



Remark 5 Besides, if we consider the 2 -dimensional Navier-Stokes equations with the repulsive 
force in radial symmetry with 9 — f/2, 



f 

Pt + Upr + pUr + -pu = 0, 



p {Ut + UUr) — -\ ^ /J" p(i, s)sds + [K,p^/'^\Ur + {np^/'^){Urr + "Wr — 



the special solutions are: 



p{t,r) = 



1 



(T - Ct)^y 



T-Ct 



j/(z)+iy(z) + = 0, y(0) = a ^ 0, j)(0) = 0, 



(33) 



(34) 
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